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^ ■ Abstract 



The paper contains the Weyl formula for the counting function of the interior 
transmission problem when the latter is parameter-elliptic. Branching billiard tra- 
jectories are constructed, and the second term of the Weyl asymptotics is estimated 
from above under some conditions on the set of periodic billiard trajectories. 
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a^■ 1 Interior transmission eigenvalues, 
oo 



Let O & he an open bounded domain with C°° boundary dO. The classical Weyl 
formula 



for the counting function A^(A) (number of eigenvalues whose absolute values do not exceed 
^ I A) is well known for the Dirichlet or Neumann Laplacian in O. Here V{0) is the volume 

I of the domain and Ud is the volume of the unit ball in R'^. This paper concerns the Weyl 

formula for the interior transmission eigenvalues (ITE) which were introduced by A.Kirsch 
in [T7] in connection with an inverse scattering problem for the reduced wave equation 
and further studied by D.Colton and P.Monk |^. The anisotropic interior transmission 
problem (discussed below) was introduced in [8|, see the review [9] for more references. 
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Let us recall the definition of ITE. The values of A G C for which the following problem 



-Au-Xu = 0, xeO, ueH\0), 

-VAVv - Xn{x)v = 0, xeO, v e H^{0), 



(2) 



u — V = 0, X e do, 

|H _ |3L = 0, xedO. ^ ' 

has a non-trivial solution are called the interior transmission eigenvalues. Here H'^{0) 
is the Sobolev space, A{x), x G (9, is a smooth symmetric elliptic {A = > 0) matrix 
with real-valued entries, n{x) is a smooth function, u is the outward normal to dO, and 
the CO- normal derivative is defined as follows 

V = V ■ A\v. 



du 

We will mostly be concerned with the cases d = 2,3, but all the results below can be 
automatically carried over to any dimension d. 

There are many papers on the Weyl formula for general elliptic boundary value prob- 
lems (see review [1] for references). The particular feature of the problem under consider- 
ation is that it is not symmetric. The formally conjugate problem has different boundary 
conditions: 

u + V = 0, X e do, . 
^ + — = QxGdO 

du ^ duA ' ^ 

The spectrum of the problem ([2]),(j3]) is not always discrete. Examples when the eigenvalues 
fill the whole complex plane can be found in p.8j. It was also shown there that ITE form 
a discrete set when the problem is parameter-elliptic. Conditions for parameter-ellipticity 
were also described in [12], and they will be formulated below in the case of real A and n. 
This paper contains a justification of the Weyl formula for ITE when the problem (12]),(j3]) 
is parameter-elliptic. In particular, this implies that the set of ITE is infinite for general 
parameter-elliptic problems. The infiniteness of the set of ITE was known in many cases, 
see [4J, [5j and references there. When A = Id, some estimates on A^(A) for ITE can be 
found in [21j,[T5. Before we proceed with the main result on the Weyl formula, we would 
like to show that ITE play the same role for the transmission scattering problem as the 
eigenvalues of the Dirichlet or Neumann Laplacian play for the scattering by an obstacle 
with the corresponding (Dirichlet or Neumann) boundary conditions. 

A connection between the counting function A^(A) for the Dirichlet or Neumann Lapla- 
cian and the total scattering phase (which is defined as — aj:gdetS{k), where S{k) is the 
scattering matrix) was established in 1978, see |19] . |16] . It was shown that 

-— argdet5'(A:) = A^(A)(1 + 0(A-")), e > 0, \ = k^^oo. 

This formula suggests that A^(A) can perhaps be interpreted as the rotation number of 
the scattering phase (number of rotations of exp{i arg detS{k)) on the unit circle when 
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k — !■ oo). At present, there is a deeper understanding of the relation between A^(A) and 
the rotation number. 

Recall that the operator S{k), A; > 0, is unitary and its eigenvalues z = Zj{k) belong 
to the unit circle. Consider two conditions. 

A. —k^ is an eigenvalue of the Dirichlet or Neumann Laplacian. 

B. 2 = 1 is an eigenvalue of S-matrix {S = S{k)) for the scattering problem with the 
same boundary condition and A; > 0. 

It is easy to see that B implies A. Indeed, the unitary matrix S{k) has the form S{k) = 
I+^F, where F is the integral operator on L^(5''^"^) (S**^"^ is the unit sphere) whose kernel 
is the scattering amplitude. If Sfi{9) = fj.{0), 9 G S'^^^ , then F/u = 0. Let i(){k,u,x) = 
^ik{ui,x) _|_ ^^^^ ^ ^ S'^^^ , be the solution of the scattering problem. Consider function 
u{k^ x) = Jgd-i il^^{uj)dSn. Function u satisfies the homogeneous boundary condition since 
ip satisfies it. Since F/i = 0, the outgoing wave Jgd-i '^^J'{i^)dSf^ has zero amplitude, and 
therefore it is equal to zero identically. Thus u{k,x) = fgd-i e'^''^'^'^'^ fi{u)dS^. Obviously, 
this function satisfies the Helmholtz equation in the whole space, i.e., u is an eigenfunction 
of the interior problem with the eigenvalue —k"^. 

This implication from B to A was first noted in papers on inverse scattering problem 
(see [S]) were slightly different terminology was used: injectivity of the far field operator 
F implies A. In terms of S'-matrix, it can be found in [TU]. 

The inverse implication (from A to B) holds in the case of a ball, but not for general 
domains, see However, for general domains O, B implies A in a weaker sense: 

— A;q is an eigenvalue of the Dirichlet Laplacian if and only if there exists an analytic in 
k G {ko — 6, ko), 6 > 0, eigenvalue z = Zi{k) of the S'-matrix such that 

lim arg2;j(A;) = 2n + 0. (5) 

In connection with the last relation, let us note that the eigenvalues z = {zj{k)} of 
the S'-matrix S{k) form an analytic manifold while z 1. Point z = 1 is an essential 
point of the spectrum of the operator S{k) (limiting point for the set of eigenvalues) and 
the structure of the manifold in a neighborhood of this point is much more complicated. 
In particular, it may happen that an eigenvalue approaches z = 1 as k k^, but z = 1 
is not an eigenvalue of S'(A;o). It is also essential that we have — > A;o — in ([5]), but not 
k ko + 0. 

The study of the connection between properties A and B was initiated by E. Doron, 
U. Smilansky, [TU]. Relation was justified by J. Eckmann and C. Pillett in [TT] in 
the case of the Dirichlet boundary condition. Paper |12| contains a similar result for the 
Neumann boundary condition. 

All the relations between the Dirichlet /Neumann problem and scattering problem men- 
tioned above can be easily reformulated for ITE and the (exterior) transmission scattering 
problem. The latter problem is stated as follows. 

-Au-Xu = 0, xeR\0, u = e'''^''''^'^ + 'ilJsc{x,k,u), 

-VAVv - Xn{x)v = 0, X G C, ^ ^ 
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u — V = 0, X G do, 
where A = fc^, ipse satisfies the radiation conditions: 

Proposition 1.1. If z = 1 is an eigenvalue of the scattering matrix Str{k), k > 0, for 
the transmission scattering problem then k"^ is one of the ITE. 

The proof of this statement is the same as in the case of the Dirichlet or Neumann 
boundary conditions. The inverse statement is also vahd for the transmission problem 
when is a ball, A = a Id, a and n are constant. We believe that the result on the weak 
implication from B to A also holds for the transmission problem (for arbitrary domains 
O), but it has not been proved yet. 

Let us recall conditions on A, n which guarantee the parameter-ellipticity of interior 
transmission problem ([2]), In this paper we assume that A and n are real- valued. Let 
us fix an arbitrary point x^ G dO and choose a new orthonormal basis {ej}, 1 < j < d, 
centered at the point x^ with = u, where u is the normal to the boundary at the point 
x". The vectors ei, ed-i belong to the tangent plane to dO at the point xq. Let y be 
the local coordinates defined by the basis {ej}, and let C = C{x^) be the transfer matrix, 
i.e., y = C{x — a;°). 

We fix the point x = x° in equations ([2]), ([3]) and rewrite the problem in the local 
coordinates y. Then we get the following problem with constant coefficients in the half 
space yd > : 

-AyU^- k'^u = 0, yd> 0, 

-VyAVyV - Pn{x°)v = 0, > 0, ^ ' 

u-v = 0, yd = 0, 

£-li = 0' ^'^ = 0- 

Here 

A = A{x'^) = CA{x'^)C*. 

The entries of the matrix A = ( equal to aij = ej ■ A{x^)ei. The co-normal 

derivative in the boundary condition equals ■ AVy. 
The following result can be extracted from |18) . 



Theorem 1.2. Let A{x) > 0,n{x) > for x E O and the following conditions hold for 
allx^ G dO: ifd = 2, then 

a2,2ri(x°) -1^0, and detA(x°) ^1, x° G dO; 

if d = 3, then a3^3n{x^) — 1 7^ and 

det ( ~ ^""'''^^ ~ ^ ~ A>0, x" G dO. 

\ a3,3«2,i - 01,302,3 03,302,2 - (02,3) - 1 / 
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Then the interior transmission problem ^ is parameter- elliptic for A in any 
sector of the complex X-plane which does not contain either of the rays and Its 
eigenvalues form a discrete set. 

Moreover, if a = 1, then the problem is parameter- elliptic also on Here 

a = sgn[(a2,2n(x°) - l)(detA(x°) - 1)], x° G dO, 

if d = 2. If d = 3, then a = 1 when the matrix 

a3,3n{x°) - 1 \ 

03,301,1 - (01,3)^ - 1 03,301,2 - 01,302,3 , x° e do, 

03,302,1 - 01,302,3 03,302,2 - (02,3)^ ^ 1 / 

is sign- definite and a = —1 otherwise. 



2 Weyl asymptotics and branching billiards. 

Let us split the set {A}, A G C, of ITE into two subsets {A+} and {A^} where ReA+ > 0, 
ReA~ < 0. We enumerate the ITE A^,A~ in increasing order of |A„| and denote by 
N'^{t),N^{t) the counting functions for {A^},{A^}. Similarly, we enumerate the whole 
set of the ITE in increasing order of |A„| and denote by N{t) the counting functions for 
{A„}. 

Theorem 2.1. Let assumptions of Theorem \1.2\ hold (with arbitrary a = ±1). Then 

1 ) There is at most a finite number of ITE inside any sector of complex X-plane which 
does not contain either of the rays and R_ . 

2) The Weyl formula holds for {A^}; 

N+{t) ~ at'^/^ t 00, (10) 

where 

(2'kYJo\ idetA{x)y/y ^ ' 

In particular, if A = aid, where a > 0, and n is also constant, then 

3) If a = 1 then the set {A^} is finite. If a = —1 then there exists M > such that 

N-{t)<Mt, t^oo, d = 2; N"{t) < Mtlnt, t ^ 00, d = 3. 
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Remark If n{x) = a, A{x) = aid, where a > 0, a ^ 1, then the substitution u — v = 
Ml, au — V = Vi reduces the problem ([2]), ([3]) to the Dirichlet problem for Ui and the 
Neumann problem for Vi. Hence in this case, the set {A„} coincides with the union of the 
Dirichlet and Neumann Laplacians. 

Proof. The first statement of the theorem and the first part of the last statement 
(concerning a = 1) are immediate consequences of the parameter-ellipticity of the problem 
established in Theorem 11.21 The second statement follows from Theorem 11.21 and results 
obtained in [2j,[3j where the main term of the Weyl asymptotics is justified for non- 
symmetric parameter-elliptic system if the eigenvalues of the main symbol A of the system 
belong to Formula (fTOl) holds in this case with 

[2Tlf Jo Jjid 

where N{l,x,^) is the number of eigenvalues of A whose absolute values do not exceed 
one. In our case, 

e 

and 



A 



which implies ( fTTl) . 

The part of the third statement concerning the case a = —1 is a consequence of 
Theorem 2 from [2],|3]. The latter theorem estimates the number of the eigenvalues of 
the problem in a sector, which does not contain the eigenvalues of the principle symbol 
of the operator. □ 

It is well known that the remainder term in the classical Weyl asymptotics depends 
on the set of the closed billiard trajectories (periodic trajectories of the corresponding 
Hamiltonian system). The second term of the asymptotics was justified by V. Ivrii [15] 
for boundary value problems for the Laplacian under the condition that the set of periodic 
trajectories has zero measure. D. Vasiliev |22] and Yu. Safarov [20] extended this result 
to the case of self-adjoint elliptic system where the billiard trajectories are not defined 
uniquely after a ray hits the boundary (see figure [T^). We can not use the latter results 
directly since the problem ([2]), ([3]) is not symmetric, but we will combine these results with 
the theory of s-numbers to obtain an estimate for the second term of the asymptotics from 
above. Let us describe the branching trajectories for our problem. 

Consider two Hamiltonian flows G*, z = 1,2, with Hamiltonians 



/ii(x,0 = iei, h2{x,0 = xi^,^'A{x)i. xeO, ^eR' {=r{0)). (12) 

We assume that hi{x,^) ^ h2{x,^),\/x E O,^ E R'^\{0}, which is equivalent to the 
following: 
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Assumption 2.2. Matrix j^^A — Id is positive-definite or negative- definite when x & O. 

When a ray which corresponds to one of these Hamiltonians comes to the boundary 
do it creates two reflected rays. One of the reflection angles corresponds to the "billiard 
law": reflection angle is equal to the incident angle. This ray corresponds to the same 
Hamiltonian as the Hamiltonian of the incident ray. The second reflection ray is deflned 
buy Snell's law. One can determine both reflected rays as follows. If the Hamiltonian of 
the incident ray is equal to one (it is constant along the trajectory), then the directions 
of the reflected rays are deflned by the relations /ii = /i2 = 1 at the point of reflection 
(and it may happen that one of these equations does not have a solution in which case 
there is only one reflected ray). 

One needs to work with branching trajectories in order to construct quasi-modes 
(approximate solutions of the problem). A branching trajectory splits every time when it 
has a chance to split after the reflection from the boundary. After the initial Hamiltonian 
hifj (to start the trajectory) is chosen, the branching trajectory is defined uniquely by 
its initial data {y,ri), y & O, 77 G R^, hi^ly,!]) = 1. The initial data form a {2d — 1)- 
dimensional manifold. For our purpose (an estimate on the second term in the Weyl 
asymptotics), we need another object, billiard trajectories. These trajectories (x*,^*), t > 
0, do not split. The point moves according to one of the Hamiltonian flows and the 
Hamiltonian can be changed after each reflection from the boundary. Thus each initial 
data usually defines infinitely many billiard trajectories. 

A billiard trajectory is called a dead-end trajectory if the ray touches the boundary 
or there are infinitely many refiections on a finite time interval. 

Assumption 2.3. The measure of the dead-end trajectories (i.e., the {2d — 1)-Lebesgue 
measure of the set of initial data for the dead-end trajectories) is zero. 

Consider a periodic billiard trajectory with a period T and initial data (?/0)^o)- This 
trajectory is called absolutely periodic if for each {y,ri) in a p-neighborhood of (yo)^o) a 
billiard trajectory (x*,^*), t>0, with the initial data {y,T]) can be chosen in such a way 
(by repeating the same pattern of refiections) that the trajectory at time t = T is located 
in 0{p°°), p — > 0, neighborhood of the starting point {yji]). The latter means that the 
function (x*,^*) — {y,ri) has zero of infinite order at the point {yojfjo). The initial data 
{yo,Vo) is called absolutely periodic if at least one absolutely periodic trajectory has this 
data. 

Assumption 2.4. {2d — 1)-Lebesgue measure of absolutely periodic data is zero. 
Theorem 2.5. Let the assumptions of Theorem \1.2\ and Assumptions\2.2^2.4\ hold. Then 



there is a constant C such that 

N{t) < af^/^ + Ct(^~i)/2 ^ o(t('^-i)/2), t 00, (13) 
where a is defined in fTl\) . 
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Figure 1: Strongly periodic branching trajectory of the interior problem and the corre- 
sponding trajectory of the exterior problem. 



Remark. Assumption 12.41 can be weakened, see 
Proof. Denote by the weighted space LF'{0) x Lp'iO) with the weight which 
corresponds to the following scalar product: 

(■Ui, f i) ■ (■U25 "^2) = / UiU2dx + I n{x)viV2dx. (14) 



Consider the following operator L in which corresponds to problem ([2]), (jS]): the 
mapping L is defined by 

L{u, v) = {-Au, -^VA{x)Vv), (15) 
n[x) 

and the domain consists of vectors {u,v), u,v & H'^{0), such that u,v satisfy (|7]). 
Operator L*, adjoint to L with respect to the scalar product (HM . can be obtained from 
L if we replace the minus signs in both relations ([71) by the plus signs. 

Let us fix a real 7 which is not an ITE, i.e., the operators (L — 7/) and (L* — 7/) are 
invertible. Consider the operator M = {L — ^I)*{L — 7J) defined on {u, v),u,v G iJ^((9), 
such that {u,v) G Dl, {L — '-f){u,v) G D^*. Operator M in with the domain Dm 
described above is self-adjoint. One can rewrite the equation Mz = Xz, z = {u, v), in the 
form of the boundary value problem: 

{-A-j)^u = Xu, xeO, 
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u — V = 0, X e do, 
|H _ |!L = 0, xedO, 

(A + 7Tw+(^VAV + 7)^^ = 0, xedO, (17) 

au Ova 

If the assumptions of Theorem 11.21 hold, then problem f lT6|) . (fT7|) is elliptic. Indeed, one 
can check that the Shapiro-Lopatinskii condition for problem (fT6|) . (fT7|) coincides with 
the ones for the problem ([2]), ([3]). Another option is to note that the standard elliptic 
a priory estimates are valid for the solutions of inhomogeneous problem ( fT6l) . ( fTT|) with 
inhomogeneities in ( IT6l) (since they are valid for operators (L — 7/) and (L* — 7/)) and 
the latter implies the ellipticity of the boundary value problem. Since problem (fT6|) . 
( IT7|) is symmetric, we can apply the results of D. Vassiliev [22| Th.1.2] and Yu. Safarov 
||20| Th.1.1] which provide an estimate on the second term of the Weyl asymptotics for 
symmetric elliptic systems. Namely, it follows that if the assumptions of Theorem II. 21 and 
Assumptions 12. 2112.41 hold, then 

N{t) = ati + Ct'^ + o{t^), t^oo, (18) 

where N{t) is the counting function for operator (L — ^I)*{L — 7/) and a is defined in 
( ITT]) . Operators [(L — jiy^L — 7/)]^^ and (L — 7/)^^ are compact, and therefore (e.g., 
see [131 Lemma 3.3]) 

|A„((L - 7/)-^)P < A„([(L - 7/)*(^ - ll)]-'), (19) 

where the eigenvalues A„ are enumerated in the order of the decay of |A„| (contrary to the 
eigenvalues of operators M or L for which |A„| — 00 as n — )■ 00). Inequality fll9p implies 
that 

|A„((L-7/))r>A„([(L-7/)*(L-7/)]), 
and this together with f]T8]) justify the statement of the theorem. □ 



3 Strongly periodic branching trajectories. 

Consider again the branching billiard trajectories of the ITE problem. In the previous 
section, we considered periodic billiard trajectories where the point moves (after hitting 
the boundary) along one of the possible reflected rays. Now we assume that the trajectory 
splits every time when it has a chance to split. We will call such a trajectory strongly 
periodic branching billiard trajectory if it consists of a finite number of segments (see 
Figure [1^). Under some additional condition on the optical lengths of some parts of the 
trajectory, one can use these trajectories to construct quasi-modes for the ITE problem, 
i.e., to construct functions which almost (as k = kn ^ 00) satisfy the equations and 
boundary conditions of the ITE problem. 



9 



Figure 2: Periodic branching trajectories of the exterior transmission problem which are 
not related to the periodic branching trajectories of the interior problem. There is only 
one incident ray in a), and three incident rays in b). 



One also can construct quasi-modes for the exterior transmission problem. The corre- 
sponding billiard trajectories are defined by the Hamiltonian flow outside of the obstacle 
with the Hamiltonian hi{x,^) = \^\ and the Hamiltonian flow inside of the obstacle with 



the Hamiltonian h2{x, ^) = y After hitting the boundary, each ray splits into 

a reflected and refracted parts. Such a trajectory is called periodic if it consists of a finite 
number of segments (exterior rays can be semi-infinite, see Figures [Tb and [2]). 

While the periodic trajectories of the exterior problem may not be related to strongly 
periodic trajectories of the interior problem (Figure [2]), each strongly periodic branching 
trajectory of the interior problem always define a corresponding trajectory of the exterior 
problem. In order to define such exterior trajectory, one needs to replace each segment 
of the interior trajectory which corresponds to the Hamiltonian hi{x,^) = \^\ by its 
complement on the line containing this segment (see Figure [T]) . One also can extend to 
Bf'\0 the parts of the interior quasi-mode which correspond to hi, and then we omit 
the interior parts of the quasi-mode related to hi. This will provide a quasi-mode of 
the exterior problem which has the following property. Its exterior part can be smoothly 
extended inside of the obstacle as a solution of the Helmholtz equation in the whole 
space (where the extension is the omitted part of the interior mode). The latter is an 
important feature in the relation between the ITE and the eigenvalues of the S'-matrix. 
In particular, this feature provides an alternative proof that property B implies property 
A which is discussed in the first section. 

Acknowledgment. The authors are very grateful to F. Cakoni, D. Colton, B. Gutkin, 
H. Haddar, V. Ivrii, A. Kirsch, Yu. Safarov, D. Vassiliev for useful discussions. 
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